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ABSTRACT 


A  14.25-percent-thick,  natural-laminar-flow  airfoil,  the  S406,  for  rotorcraft  applica¬ 
tions  has  been  designed  and  analyzed  theoretically  and  verified  experimentally  in  The  Penn¬ 
sylvania  State  University  Low-Speed,  Low-Turbulence  Wind  Tunnel.  The  two  primary 
objectives  of  high  maximum  lift  and  low  profile  drag  have  been  achieved.  The  constraint  on 
the  airfoil  thickness  has  been  satisfied,  but  the  one  on  the  pitching  moment  has  not.  The  air¬ 
foil  exhibits  a  docile  stall.  Comparisons  of  the  theoretical  and  experimental  results  generally 
show  good  agreement. 


INTRODUCTION 


Almost  all  airfoils  in  use  on  rotorcraft  today  were  developed  under  the  assumption  that 
extensive  laminar  flow  is  not  likely  on  a  rotor.  (See  ref.  1,  for  example.)  For  the  present 
application,  however,  given  the  relatively  low  Reynolds  numbers  and  the  precision  blade 
manufacturing  technique  being  employed,  the  achievement  of  laminar  flow  warrants  explora¬ 
tion. 


The  airfoil  designed  under  the  present  effort  is  intended  for  the  rotor  of  a  small  heli¬ 
copter.  To  complement  the  design  effort,  an  investigation  was  conducted  in  The  Pennsylvania 
State  University  Low-Speed,  Low-Turbulence  Wind  Tunnel  (ref.  2)  to  obtain  the  basic,  low- 
speed,  two-dimensional  aerodynamic  characteristics  of  the  airfoil.  The  results  have  been 
compared  with  predictions  from  the  method  of  references  3  and  4  and  from  the  method  of  ref¬ 
erence  5. 


SYMBOLS 

Values  are  given  in  both  SI  and  U.S.  Customary  Units.  Measurements  and  calcula¬ 
tions  were  made  in  U.S.  Customary  Units. 

Cn  pressure  coefficient,  — — — 

Too 

c  airfoil  chord,  mm  (in.) 

cc  section  chord-force  coefficient, 

cd  section  profile-drag  coefficient, 

cnsina  +  cccosa 

cd’  point  drag  coefficient  (ref.  6) 

c i  section  lift  coefficient,  cn/ cosa  -  cdtana 


#cpd(D 

J  cd'd^)  ,  except  post  stall, 

Wake 


h 

M 

P 

q 

R 

s 

t 

X 

y 

z 

a 

Subscripts: 


section  pitching-moment  coefficient  about  quarter-chord  point, 

-icp(ra25)d©+K@d© 

section  normal-force  coefficient,  -  j)Cpd(^ 
horizontal  width  in  wake  profile,  mm  (in.) 
free-stream  Mach  number 

'y 

static  pressure,  Pa  (lbf/ft  ) 

'•y 

dynamic  pressure,  Pa  (lbf/ft  ) 

Reynolds  number  based  on  free-stream  conditions  and  airfoil  chord 

arc  length  along  airfoil  surface,  mm  (in.) 

airfoil  thickness,  mm  (in.) 

airfoil  abscissa,  mm  (in.) 

model  span  station,  y  =  0  at  midspan,  mm  (in.) 

airfoil  ordinate,  mm  (in.) 

angle  of  attack  relative  to  x-axis,  deg 


/  local  point  on  airfoil 

11  lower  limit  of  low-drag  range 

max  maximum 

min  minimum 

S  separation 

T  transition 

ul  upper  limit  of  low-drag  range 

0  zero  lift 
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oo 


free-stream  conditions 


Abbreviations: 

L.  lower  surface 

NACA  National  Advisory  Committee  for  Aeronautics 

S.  boundary-layer  separation  location,  xs/c 

T.  boundary-layer  transition  location,  xT/c 

U.  upper  surface 


AIRFOIL  DESIGN 
OBJECTIVES  AND  CONSTRAINTS 

The  airfoil  design  specifications  are  contained  in  table  I.  Two  primary  objectives  are 
evident.  The  first  objective  is  to  achieve  a  maximum  lift  coefficient  of  1 .30  at  a  Mach  number 
of  0.30  and  a  Reynolds  number  of  1.14  x  106.  A  requirement  related  to  this  objective  is  that 
the  maximum  lift  coefficient  not  decrease  significantly  with  transition  fixed  near  the  leading 
edge  on  both  surfaces.  In  addition,  the  airfoil  should  exhibit  docile  stall  characteristics.  The 
second  objective  is  to  obtain  low  profile-drag  coefficients  from  a  lift  coefficient  of  0.20  at  a 
Mach  number  of  0.59  and  a  Reynolds  number  of  2.12  x  106  to  a  lift  coefficient  of  0.70  at  a 
Mach  number  of  0.46  and  a  Reynolds  number  of  1.63  x  106. 

Two  major  constraints  were  placed  on  the  design  of  the  airfoil.  First,  the  zero-lift 
pitching-moment  coefficient  must  be  no  more  negative  than  -0.05.  Second,  the  airfoil  thick¬ 
ness  must  equal  14.25-percent  chord. 


3 


PHILOSOPHY 


Given  the  above  objectives  and  constraints,  certain  characteristics  of  the  design  are 
apparent.  The  following  sketch  illustrates  a  drag  polar  that  meets  the  goals  for  this  design. 


Sketch  1 


The  desired  airfoil  shape  can  be  traced  to  the  pressure  distributions  that  occur  at  the  various 
points  in  sketch  1 .  Point  A  is  the  lower  limit  of  the  low-drag  range  of  lift  coefficients;  point  B, 
the  upper  limit.  The  profile-drag  coefficient  at  point  B  is  not  as  low  as  at  point  A,  unlike  the 
polars  of  many  laminar-flow  airfoils  where  the  drag  coefficient  within  the  laminar  bucket  is 
nearly  constant.  (See,  for  example,  ref.  7.)  This  characteristic  is  related  to  the  elimination  of 
significant  (drag-producing)  laminar  separation  bubbles  on  the  upper  surface  for  the  design 
range  of  Reynolds  numbers.  (See  ref.  8.)  The  drag  coefficient  increases  rapidly  outside  the 
low-drag,  lift-coefficient  range  because  boundary-layer  transition  moves  quickly  toward  the 
leading  edge  with  increasing  (or  decreasing)  lift  coefficient.  This  feature  results  in  a  leading 
edge  that  produces  a  suction  peak  at  higher  lift  coefficients,  which  ensures  that  transition  on 
the  upper  surface  will  occur  very  near  the  leading  edge.  Thus,  the  maximum  lift  coefficient, 
point  C,  occurs  with  turbulent  flow  along  the  entire  upper  surface  and,  therefore,  should  be 
relatively  insensitive  to  roughness  at  the  leading  edge. 
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From  the  preceding  discussion,  the  pressure  distributions  along  the  polar  can  be 
deduced.  The  pressure  distribution  at  point  A  should  look  something  like  sketch  2. 


Sketch  2 


To  achieve  low  drag,  a  favorable  pressure  gradient  is  desirable  along  the  upper  surface  to 
about  30-percent  chord.  Aft  of  this  point,  a  short  region  having  a  shallow,  adverse  pressure 
gradient  (“transition  ramp”)  promotes  the  efficient  transition  from  laminar  to  turbulent  flow 
(ref.  9).  The  transition  ramp  is  followed  by  a  slightly  convex  pressure  recovery.  The  specific 
pressure  recovery  employed  represents  a  compromise  between  maximum  lift,  drag,  pitching 
moment,  stall  characteristics,  and  drag  divergence.  The  steep,  adverse  pressure  gradient  aft  of 
about  90-percent  chord  is  a  “separation  ramp,”  originally  proposed  by  F.  X.  Wortmann,1 
which  confines  turbulent  separation  to  a  small  region  near  the  trailing  edge.  By  constraining 
the  movement  of  the  separation  point  at  high  angles  of  attack,  higher  lift  coefficients  can  be 
achieved  with  little  drag  penalty.  This  feature  has  the  added  benefit  of  promoting  docile  stall 
characteristics.  (See  ref.  10.) 

Along  the  lower  surface,  the  pressure  gradient  is  initially  adverse  and  then  approxi¬ 
mately  zero  to  about  65-percent  chord.  Thus,  transition  is  imminent  over  the  entire  forward 
portion  of  the  lower  surface.  (See  ref.  11.)  This  concept  allows  a  wide  low-drag  range  to  be 
achieved  and  increases  the  loading  in  the  leading-edge  region.  The  forward  loading  serves  to 
balance,  with  respect  to  the  pitching-moment  constraint,  the  aft  loading,  both  of  which  con¬ 
tribute  to  the  achievement  of  the  specified  maximum  lift  coefficient  and  low  profile-drag  coef¬ 
ficients.  This  region  is  followed  by  a  transition  ramp  and  then  a  roughly  linear  pressure 
recovery.  The  pressure  recovery  must  begin  farther  forward  than  optimum  for  low  drag  and 
the  constrained  pitching  moment  to  alleviate  separation  at  lower  lift  coefficients,  especially 
with  transition  fixed  near  the  leading  edge. 


'Director,  Institute  for  Aerodynamics  and  Gas  Dynamics,  University  of  Stuttgart,  Germany,  1974-1985. 
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The  amounts  of  pressure  recovery  on  the  upper  and  lower  surfaces  are  detennined  by 
the  airfoil-thickness  and  pitching-moment  constraints. 

At  point  B,  the  pressure  distribution  should  look  like  sketch  3. 


Sketch  3 


No  suction  peak  exists  at  the  leading  edge.  Instead,  a  rounded  peak  occurs  aft  of  the  leading 
edge,  which  allows  some  laminar  flow  although  not  to  the  extent  of  point  A. 


EXECUTION 

Given  the  pressure  distributions  previously  discussed,  the  design  of  the  airfoil  is 
reduced  to  the  inverse  problem  of  transforming  the  pressure  distributions  into  an  airfoil  shape. 
The  Eppler  Airfoil  Design  and  Analysis  Code  (refs.  3  and  4)  was  used  because  of  its  unique 
capability  for  multipoint  design  and  because  of  confidence  gained  during  the  design,  analysis, 
and  experimental  verification  of  many  other  airfoils.  (See  ref.  12,  for  example.) 

The  airfoil  is  designated  the  S406.  The  airfoil  shape  and  coordinates  are  available 
from  Airfoils,  Incorporated.  The  airfoil  thickness  is  14.25-percent  chord,  which  satisfies  the 
design  constraint. 
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THEORETICAL  PROCEDURE 


The  theoretical  results  are  predicted  using  the  method  of  references  3  and  4 
(PROFIL07),  commonly  known  as  the  Eppler  code,  and  the  method  of  reference  5 
(MSES  3.0).  Critical  amplification  factors  of  11  and  9  were  specified  for  the  computations 
using  the  method  of  references  3  and  4  and  the  method  of  reference  5,  respectively.  Because 
the  maximum  lift  coefficient  computed  by  the  method  of  references  3  and  4  is  not  always  real¬ 
istic,  an  empirical  criterion  has  been  applied  to  the  computed  results.  The  criterion  assumes 
the  maximum  lift  coefficient  has  been  reached  if  the  drag  coefficient  of  the  upper  surface 
reaches  a  certain  value  that  is  a  function  of  the  Reynolds  number  and  the  wind-tunnel  facility. 

Because  the  free-stream  Mach  number  for  all  wind-tunnel  test  conditions  did  not 
exceed  0.2,  the  flow  can  be  considered  incompressible  for  the  purpose  of  comparing  the  theo¬ 
retical  and  experimental  results.  This  allows  the  (incompressible)  conformal-mapping 
(design)  method  of  references  3  and  4  and  the  fast,  subcritical  flow  solver  of  reference  5  to  be 
used. 


EXPERIMENTAL  PROCEDURE 
WIND  TUNNEL 

The  Pennsylvania  State  University  Low-Speed,  Low-Turbulence  Wind  Tunnel  (ref.  2) 
is  a  closed-throat,  single-return,  atmospheric  tunnel  (fig.  1).  The  test  section  is  101.3  cm 
(39.9  in.)  high  by  147.6  cm  (58.1  in.)  wide  (fig.  2).  Electrically  actuated  turntables  provide 
positioning  and  attachment  for  the  two-dimensional  model.  The  turntables  are  flush  with  the 
top  and  bottom  tunnel  walls  and  rotate  with  the  model.  The  axis  of  rotation  coincided  approx¬ 
imately  with  the  midchord  of  the  model,  which  was  mounted  vertically  between  the  turnta¬ 
bles.  The  gaps  between  the  model  and  the  turntables  were  sealed.  The  turbulence  intensity  in 
the  test  section  is  approximately  0.05  percent  at  46  m/s  (150  ft/s). 


MODEL 

The  aluminum,  wind-tunnel  model  was  fabricated  by  Advanced  Technologies,  Incor¬ 
porated,  Newport  News,  Virginia,  using  a  numerically  controlled  milling  machine.  The  model 
had  a  chord  of  457.04  mm  (17.994  in.)  and  a  span  of  107.95  cm  (42.50  in.)  and,  thus,  extended 
through  both  turntables.  Upper-  and  lower-surface  orifices  were  located  to  one  side  of  mid¬ 
span  at  the  staggered  positions  listed  in  table  II.  All  the  orifices  were  0.51  mm  (0.020  in.)  in 
diameter,  except  the  trailing-edge  orifice,  which  was  0.25  mm  (0.010  in.)  in  diameter,  with 
their  axes  perpendicular  to  the  surface.  The  surfaces  of  the  model  were  sanded  to  ensure  an 
aerodynamically  smooth  finish.  The  measured  model  contour  was  within  0.13  mm  (0.005  in.) 
of  the  prescribed  shape. 
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WAKE-SURVEY  PROBE 


A  total-  and  static-pressure,  wake-survey  probe  (fig.  3)  was  mounted  from  the  top  tun¬ 
nel  wall  (fig.  2).  The  probe  was  positioned  53  cm  (21  in.)  from  the  top  wall  (i.e.,  about  mid¬ 
span)  and  automatically  aligned  with  the  wake-centerline  streamline.  A  traverse  mechanism 
incrementally  positioned  the  probe  to  survey  the  wake.  The  increment  was  1.27  mm 
(0.050  in.)  for  traverses  less  than  254.0  mm  (10.00  in.)  and  2.54  mm  (0.100  in.)  for  longer  tra¬ 
verses,  which  were  occasionally  required  near  the  maximum  angle  of  attack.  The  tip  of  the 
probe  was  located  0.83  chord  downstream  of  the  trailing  edge  of  the  model. 


INSTRUMENTATION 

Basic  tunnel  pressures  and  the  wake  pressures  were  measured  with  precision  transduc¬ 
ers.  Measurements  of  the  pressures  on  the  model  were  made  by  an  automatic  pressure¬ 
scanning  system  utilizing  precision  transducers.  Data  were  obtained  and  recorded  by  an  elec¬ 
tronic  data-acquisition  system. 


METHODS 

The  pressures  measured  on  the  model  were  reduced  to  standard  pressure  coefficients 
and  numerically  integrated  to  obtain  section  nonnal-force  and  chord-force  coefficients  and 
section  pitching-moment  coefficients  about  the  quarter-chord  point.  Section  profile-drag 
coefficients  were  computed  from  the  wake  total  and  static  pressures  by  the  method  of  refer¬ 
ence  6.  Wake  surveys  were  not  perfonned,  however,  at  most  post-stall  angles  of  attack,  in 
which  case,  the  profile-drag  coefficients  were  computed  from  the  nonnal-  and  chord-force 
coefficients. 

Standard,  low-speed,  wind-tunnel  boundary  corrections  (ref.  13)  have  been  applied  to 
the  data.  It  should  be  noted,  however,  that  the  pressure  distributions  themselves  are  uncor¬ 
rected.  The  wake-survey-probe  total-pressure-tube  displacement  correction  (ref.  6)  has  been 
taken  into  account. 


TESTS 

The  model  was  tested  at  Reynolds  numbers  based  on  airfoil  chord  of  0.5  x  1 06, 
0.7  x  106,  1.0  x  106,  and  1.5  x  106  with  transition  free  (smooth)  and  with  transition  fixed  by 
roughness  at  2-percent  chord  on  the  upper  surface  and  7-percent  chord  on  the  lower  surface. 
The  grit  roughness  was  sized  using  the  method  of  reference  14  and  sparsely  distributed  along 
3-mm  (0.1-in.)  wide  strips  applied  to  the  model  with  lacquer.  (See  table  111(a).)  The  model 
was  also  tested  with  a  roughness  equivalent  to  NACA  standard  roughness  (ref.  7),  which  con¬ 
sisted  of  grit  roughness  having  a  nominal  size  of  0.211  mm  (0.0083  in.)  applied  to  the  model 
with  lacquer  and  sparsely  distributed  from  the  leading  edge  to  an  arc  length  of  8-percent  chord 
on  the  upper  and  lower  surfaces.  (See  table  111(b).)  (The  grit  size  was  scaled  from  the  NACA 


standard -roughness  grit  size  by  the  ratio  of  the  model  chords  used  in  the  two  wind  tunnels: 
457.04  mm  (17.994  in.)  in  the  present  investigation  and  609.6  mm  (24.00  in.)  in  the  NACA 
tests.)  The  Mach  number  did  not  exceed  0.2  for  any  test  condition. 

It  should  be  noted  that  the  test  Mach  numbers  are  much  lower  than  the  operational  val¬ 
ues  of  the  intended  application. 

Starting  from  0°,  the  angle  of  attack  was  increased  to  post-stall  values.  The  angle  of 
attack  was  then  decreased  from  0°  to  below  that  for  zero  lift. 


DISCUSSION  OF  RESULTS 
THEORETICAL  RESULTS 
Pressure  Distributions 

The  inviscid  pressure  distributions  at  various  angles  of  attack  and  Mach  numbers  pre¬ 
dicted  using  the  method  of  references  3  and  4  are  shown  in  figure  4. 


Section  Characteristics 

The  section  characteristics  at  the  three  design  conditions  with  transition  free  and  fixed 
are  shown  in  figures  5  through  7.  Based  on  the  predictions,  all  the  design  objectives  and  con¬ 
straints  have  been  met,  except  for  the  zero-lift  pitching-moment  coefficient,  which  exceeds 
the  constraint. 


EXPERIMENTAL  RESULTS 
Pressure  Distributions 

The  pressure  distributions  at  various  angles  of  attack  for  a  Reynolds  number  of 
1 .00  x  106  and  a  Mach  number  of  0. 1 1  with  transition  free  are  shown  in  figure  8.  At  an  angle 
of  attack  of  -4.05°  (fig.  8(a)),  transition  probably  occurs  around  75-percent  chord  on  the 
upper  surface  and  near  the  leading  edge  on  the  lower  surface.  At  an  angle  of  attack  of -2.02° 
(fig.  8(a)),  which  corresponds  to  the  lower  limit  of  the  low-drag,  lift-coefficient  range,  a  short 
laminar  separation  bubble  is  evident  on  the  lower  surface  around  75-percent  chord.  As  the 
angle  of  attack  is  increased,  a  short  laminar  separation  bubble  becomes  more  evident  on  the 
upper  surface  and  moves  forward  whereas  the  bubble  on  the  lower  surface  moves  slowly  aft 
(figs.  8(a)-8(c)).  At  an  angle  of  attack  of  9.17°  (fig.  8(c)),  turbulent,  trailing-edge  separation 
occurs  on  the  upper  surface.  The  amount  of  separation  increases  with  increasing  angle  of 
attack  (figs.  8(c)  and  8(d)).  The  maximum  lift  coefficient  occurs  at  an  angle  of  attack  of 
13.20°  (fig.  8(d)).  As  the  angle  of  attack  is  increased  further,  the  separation  point  continues  to 
move  forward,  although  the  leading-edge  pressure  peak  does  not  fall  (fig.  8(e)). 
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Section  Characteristics 


The  section  characteristics  with  transition  free,  transition  fixed,  and  scaled,  NACA 
standard  roughness,  denoted  “rough,”  are  shown  in  figure  9  and  tabulated  in  the  appendix. 
For  a  Reynolds  number  of  1.00  x  106  and  a  Mach  number  of  0.11  with  transition  free 
(fig.  9(c)),  the  maximum  lift  coefficient  is  1.25.  For  a  Reynolds  number  of  1.50  x  106  and  a 
Mach  number  of  0. 17  with  transition  free  (fig.  9(d)),  the  zero-lift  pitching-moment  coefficient 
is  -0.066,  the  lower  limit  of  the  low-drag  range  of  lift  coefficients  is  about  0.15,  and  the  max¬ 
imum  lift-to-drag  ratio  occurs  at  a  lift  coefficient  of  about  0.86.  (Because  the  upper  limit  of 
the  low-drag  range  is  not  sharply  defined,  a  precise  value  for  the  upper  limit  cannot  be  given.) 

The  effects  of  Reynolds  number  on  the  section  characteristics  are  summarized  in  fig¬ 
ure  10.  In  general,  with  transition  free,  the  lift-curve  slope,  the  maximum  lift  coefficient,  the 
lower  limit  of  the  low-drag  range,  and  the  magnitude  of  the  pitching-moment  coefficients 
increase  with  increasing  Reynolds  number.  The  zero-lift  angle  of  attack  and  the  zero-lift 
pitching-moment  coefficient  are  relatively  unaffected  by  Reynolds  number.  The  profile-drag 
coefficients  generally  decrease  with  increasing  Reynolds  number.  The  airfoil  exhibits  docile 
stall  characteristics  that  become  less  docile  with  increasing  Reynolds  number. 

The  effect  of  fixing  transition  on  the  section  characteristics  is  shown  in  figure  9.  In 
general,  the  lift-curve  slope,  the  maximum  lift  coefficient,  and  the  magnitude  of  the  pitching- 
moment  coefficients  decrease  with  transition  fixed.  These  results  are  primarily  a  consequence 
of  the  boundary-layer  displacement  effect,  which  decambers  the  airfoil  because  the  displace¬ 
ment  thickness  is  greater  with  transition  fixed  than  with  transition  free.  In  addition,  the  maxi¬ 
mum  lift  coefficient  decreases  with  transition  fixed  because  the  roughness  induces  earlier 
trailing-edge  separation.  The  reduction  in  maximum  lift  coefficient  is  small,  averaging  4  per¬ 
cent  over  the  test  Reynolds  number  range.  The  zero-lift  angle  of  attack  and  the  zero-lift 
pitching-moment  coefficient  are  relatively  unaffected  by  fixing  transition. 

It  should  be  noted  that,  for  most  test  conditions,  the  Reynolds  number  based  on  local 
velocity  and  boundary-layer  displacement  thickness  at  the  roughness  location  is  too  low  to 
support  turbulent  flow.  (See  ref.  15.)  Accordingly,  to  force  transition,  the  roughness  must  be 
so  large  that  it  increases  the  displacement  thickness,  which  abnormally  decreases  the  lift  coef¬ 
ficient  and  the  magnitude  of  the  pitching-moment  coefficient  and  increases  the  drag  coeffi¬ 
cient.  Conversely,  at  low  lift  coefficients,  the  roughness  on  the  upper  surface,  which  is  sized 
for  higher  lift  coefficients,  is  too  small  to  force  transition,  resulting  in  incorrectly  low  drag 
coefficients. 

The  effect  of  the  scaled,  NACA  standard  roughness  on  the  section  characteristics  is 
shown  in  figure  9.  The  effect  is  more  severe  than  that  of  fixing  transition.  The  reduction  in 
maximum  lift  coefficient  is  much  larger,  averaging  1 8  percent  for  Reynolds  numbers  greater 
than  0.5  x  106.  For  a  Reynolds  number  of  0.5  x  106,  the  height  of  the  scaled,  NACA  standard 
roughness  is  less  than  the  critical  height  computed  using  the  method  of  reference  14.  (See 
table  III.)  Accordingly,  the  effects  of  this  roughness  for  a  Reynolds  number  of  0.50  x  106  are 
inconsistent  (fig.  9(a)).  It  should  be  remembered  that  the  effect  of  roughness  is  proportional  to 
the  ratio  of  the  roughness  height  to  the  boundary-layer  thickness.  Because  the  height  of  the 
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scaled,  NACA  standard  roughness  and  the  airfoil  chord  are  constant,  the  effect  of  this  rough¬ 
ness  typically  increases  with  increasing  Reynolds  number  (because  increasing  Reynolds  num¬ 
ber  results  in  decreasing  boundary-layer  thickness). 

The  variations  of  maximum  lift  coefficient  and  minimum  profile-drag  coefficient  with 
Reynolds  number  are  shown  in  figures  11  and  12,  respectively.  With  transition  free,  the  max¬ 
imum  lift  coefficient  increases  with  increasing  Reynolds  number  whereas  the  minimum 
profile-drag  coefficient  decreases,  which  are  typical  trends  for  most  airfoils.  (The  maximum 
lift  coefficient  and  minimum  drag  coefficient  for  a  Reynolds  number  of  0.50  x  106  with  scaled 
roughness  are  too  high  and  too  low,  respectively,  probably  because  the  roughness  is  too  small 
to  force  transition  immediately,  as  previously  discussed.) 


COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 

Pressure  Distributions 

The  comparison  of  the  theoretical  and  experimental  pressure  distributions  at  various 
angles  of  attack  is  shown  in  figure  13.  It  should  be  noted  that  the  pressure  distributions  pre¬ 
dicted  using  the  method  of  references  3  and  4  (PROFIL07)  are  inviscid  and  incompressible, 
whereas  the  pressure  distributions  predicted  using  the  method  of  reference  5  (MSES  3.0)  as 
well  as  the  experimental  pressure  distributions  were  obtained  for  a  Reynolds  number  of 
1.00  x  106  and  a  Mach  number  of  0.11  with  transition  free.  It  should  also  be  noted  that  the 
theoretical  lift  coefficient  from  the  method  of  references  3  and  4  is  calculated  from  the  lift- 
curve  slope  and  the  angle  of  attack  relative  to  the  zero-lift  line,  whereas  the  lift  coefficient 
from  the  method  of  reference  5  and  from  the  experiment  is  derived  from  the  integrated  pres¬ 
sure  distribution.  (See  refs.  3-6.)  Thus,  at  a  given  lift  coefficient,  the  pressure  distribution 
predicted  using  the  method  of  references  3  and  4  does  not  necessarily  have  the  same  area  as 
the  measured  pressure  distribution.  It  should  be  noticed  that  the  angle  of  attack  shown  in 
figure  13  is  the  theoretical  value  from  the  method  of  references  3  and  4,  not  the  experimental 
value.  Also,  the  lift  coefficient  shown  in  this  figure  only  is  the  uncorrected  value. 

With  respect  to  the  method  of  references  3  and  4,  at  a  lift  coefficient  of  0.13 
(fig.  13(a)),  which  corresponds  to  the  lower  limit  of  the  low-drag  range,  the  pressure  coeffi¬ 
cients  and  the  pressure  gradients  agree  well,  except  where  laminar  separation  bubbles  are 
present  and  near  the  trailing  edge.  The  latter  disparity  is  probably  the  result  of  the  boundary- 
layer  displacement  effect.  At  a  lift  coefficient  of  0.67  (fig.  13(b)),  although  the  pressure  coef¬ 
ficients  do  not  match  exactly,  the  pressure  gradients  agree  reasonably  well,  again  except 
where  bubbles  are  present  and  near  the  trailing  edge.  At  the  experimental  maximum  lift  coef¬ 
ficient  (fig.  13(c)),  the  agreement  is  poor  because  the  effect  of  the  upper-surface,  trailing-edge 
separation  on  the  pressure  distribution  is  not  modelled  in  the  method  of  references  3  and  4. 

With  respect  to  the  method  of  reference  5,  at  a  lift  coefficient  of  0.13  (fig.  13(a)),  the 
pressure  coefficients  and  the  pressure  gradients  agree  remarkably  well.  The  location  of  the 
lower-surface  laminar  separation  bubble  is  predicted  well,  but  that  of  the  upper-surface  bubble 
is  aft  of  the  measured  location.  At  a  lift  coefficient  of  0.67  (fig.  13(b)),  although  the  pressure 
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coefficients  do  not  match  exactly,  the  pressure  gradients  agree  well.  The  predicted  location  of 
the  upper-surface  bubble  is  again  aft  of  the  measured  location.  At  the  experimental  maximum 
lift  coefficient  (fig.  13(c)),  the  agreement  is  poor  because  the  extensive,  upper-surface, 
trailing-edge  separation  is  not  predicted  by  the  method  of  reference  5. 


Section  Characteristics 

The  comparison  of  the  theoretical  and  experimental  section  characteristics  with  transi¬ 
tion  free  is  shown  in  figure  14.  The  maximum  lift  coefficient  estimated  using  the  previously 
discussed  empirical  criterion  applied  to  the  predictions  from  the  method  of  references  3  and  4 
(PROFIL07)  agrees  well  with  the  measurements.  The  method  of  reference  5  (MSES  3.0)  sig¬ 
nificantly  overpredicts  the  maximum  lift  coefficient.  The  method  of  references  3  and  4  gener¬ 
ally  overpredicts  the  profile-drag  coefficients  and  the  magnitude  of  the  pitching-moment 
coefficients,  although  the  agreement  improves  with  increasing  Reynolds  number.  The  method 
of  reference  5  generally  underpredicts  the  drag  coefficients  and  overpredicts  the  magnitude  of 
the  pitching-moment  coefficients  and  the  agreement  worsens  with  increasing  Reynolds  num¬ 
ber.  Both  methods  predict  the  zero-lift  angle  of  attack  and  the  zero-lift  pitching-moment  coef¬ 
ficient,  the  lower  limit  of  the  low-drag  range,  and  the  lift-curve  slope  reasonably  well, 
although  both  underpredict  the  effect  of  the  trailing-edge  separation  on  the  lift  coefficient  at 
higher  angles  of  attack. 

The  comparison  of  the  theoretical  and  experimental  section  characteristics  with  transi¬ 
tion  fixed  is  shown  in  figure  15.  In  general,  the  predicted  characteristics  show  similar  tenden¬ 
cies  as  with  transition  free,  although  the  general  agreement  is  poorer,  particularly  with  respect 
to  the  drag  coefficients,  probably  because  of  the  abnonnalities  introduced  by  the  roughness,  as 
discussed  previously. 


CONCLUDING  REMARKS 


A  14. 25 -percent- thick,  natural-laminar-flow  airfoil,  the  S406,  intended  for  rotorcraft 
applications  has  been  designed  and  analyzed  theoretically  and  verified  experimentally  in  The 
Pennsylvania  State  University  Low-Speed,  Low-Turbulence  Wind  Tunnel.  The  two  primary 
objectives  of  a  high  maximum  lift  coefficient  and  low  profile-drag  coefficients  have  been 
achieved.  The  constraint  on  the  airfoil  thickness  has  been  satisfied,  but  the  one  on  the  zero-lift 
pitching-moment  coefficient  has  not.  The  airfoil  exhibits  a  docile  stall.  Comparisons  of  the 
theoretical  and  experimental  results  generally  show  good  agreement. 
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TABLE  I.-  AIRFOIL  DESIGN  SPECIFICATIONS 


Parameter 

Objective/ 

Constraint 

Mach 

Number 

M 

Reynolds 

Number 

R 

Priority 

Minimum  lift  coefficient 

G,min 

0.10 

0.59 

2.12  x  106 

Low 

Maximum  lift  coefficient 

G,  max 

1.30 

0.30 

1.14  x  106 

High 

Lower  limit  of  low-drag, 
lift-coefficient  range  c/jj 

0.20 

0.59 

2.12  x  106 

Medium 

Upper  limit  of  low-drag, 
lift-coefficient  range  C/  u] 

0.70 

0.46 

1.63  x  106 

High 

Zero-lift  pitching-moment 
coefficient  cm0 

>  -0.05 

0.59 

2.12  x  106 

Low 

Thickness  t/c 

0.1425 

High 

Other  requirements: 

Natural-laminar-flow  airfoil 

Maximum  lift  coefficient  c /  max  independent  of  leading-edge  roughness 

Docile  stall  characteristics 
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TABLE  II.-  MODEL  ORIFICE  LOCATIONS 


[c  =  457.04  mm  (17.994  in.)] 


Upper  Surface 


Lower  Surface 


x/c 

y,  mm  (in.) 

x/c 

y,  mm  (in.) 

0.00000 

-139.96  (-5.510) 

0.00079 

-153.15  (-6.029) 

.00031 

-138.95  (-5.470) 

.00848 

-152.14  (-5.990) 

.00440 

-137.93  (-5.430) 

.02202 

-150.86  (-5.939) 

.01257 

-136.92  (-5.390) 

.04143 

-149.08  (-5.869) 

.02566 

-135.64  (-5.340) 

.06613 

-147.06  (-5.790) 

.04235 

-134.12  (-5.280) 

.09573 

-144.52  (-5.690) 

.06311 

-132.34  (-5.210) 

.12987 

-141.72  (-5.580) 

.08804 

-130.32  (-5.131) 

.16813 

-138.42  (-5.450) 

.11600 

-128.02  (-5.040) 

.21003 

-134.87  (-5.310) 

.14785 

-125.24  (-4.931) 

.25499 

-131.06  (-5.160) 

.18296 

-122.18  (-4.810) 

.30266 

-126.99  (-5.000) 

.22104 

-118.88  (-4.680) 

.35222 

-122.93  (-4.840) 

.26184 

-115.58  (-4.550) 

.40325 

-118.61  (-4.670) 

.30481 

-112.02  (-4.410) 

.45511 

-114.29  (-4.500) 

.34998 

-108.21  (-4.260) 

.50723 

-109.97  (-4.329) 

.39684 

-104.41  (-4.110) 

.55889 

-105.65  (-4.159) 

.44477 

-100.34  (-3.950) 

.60956 

-101.32  (-3.989) 

.49359 

-96.28  (-3.790) 

.65862 

-97.27  (-3.830) 

.54278 

-92.20  (-3.630) 

.70553 

-93.46  (-3.680) 

.59187 

-88.15  (-3.471) 

.75096 

-89.65  (-3.530) 

.64040 

-84.08  (-3.310) 

.79465 

-86.09  (-3.389) 

.68785 

-80.02  (-3.150) 

.83616 

-82.54  (-3.250) 

.73377 

-76.21  (-3.000) 

.87461 

-79.23  (-3.119) 

.77756 

-81.54  (-3.210) 

.91038 

-76.18  (-2.999) 

.81914 

-86.39  (-3.401) 

.94209 

-87.36  (-3.439) 

.85703 

-90.96  (-3.581) 

.96727 

-96.26  (-3.790) 

.89108 

-95.00  (-3.740) 

.98558 

-102.59  (-4.039) 

.92369 

-98.99  (-3.897) 

.99683 

-106.28  (-4.184) 

.95101 

-102.29  (-4.027) 

1.00000 

-109.20  (-4.299) 

.97255 

-104.83  (-4.127) 

.98815 

-106.61  (-4.197) 

.99681 

-107.88  (-4.247) 
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TABLE  III.-  ROUGHNESS  LOCATIONS  AND  SIZES 


(a)  Transition  fixed 


R 

Upper  surface 

Lower  surface 

x/c 

Grit 

number 

Nominal  size, 
mm  (in.) 

x/c 

Grit 

number 

Nominal  size, 
mm  (in.) 

0.5  x  106 

0.02 

60 

0.297  (0.0117) 

0.07 

30 

0.711  (0.0280) 

0.7  x  106 

80 

0.211  (0.0083) 

36 

0.589  (0.0232) 

1.0  x  106 

100 

0.150  (0.0059) 

54 

0.351  (0.0138) 

1.5  x  106 

120 

0.124  (0.0049) 

70 

0.249  (0.0098) 

(b)  Scaled,  NACA  standard  roughness 


Upper  surface 

Lower  surface 

R 

s/c 

Grit 

Nominal  size, 

s/c 

Grit 

Nominal  size, 

number 

mm  (in.) 

number 

mm  (in.) 

0.5  x  106 

0.7  x  106 

0  to 

80 

0.211  (0.0083) 

0  to 

80 

0.211  (0.0083) 

1.0  x  106 

0.08 

0.08 

1.5  x  106 
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8  m 
(25  ft) 


224-kW 

(300-hp) 

motor 


Figure  1.-  The  Pennsylvania  State  University  Low-Speed,  Low-Turbulence  Wind  Tunnel. 


Figure  2  -  Typical  airfoil  model  and  wake-survey  probe  mounted  in  test  section. 


Static-pressure  connection 
Total-pressure  connection 


Figure  3.-  Wake-survey  probe. 


(b)  M  =  0.46. 
Figure  4.-  Continued. 
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Sonic  limit 


(c)  M  =  0.59. 
Figure  4.-  Concluded. 
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PR0FIL07 
MSES  3.0 


(a)  Transition  free. 

Figure  5.-  Theoretical  section  characteristics  at  M  =  0.30  and  R  =  1.14  x  106. 


PR0FIL07 


(b)  Transition  fixed. 


K> 


Figure  5.-  Concluded. 


PR0FIL07 
MSES  3.0 


(a)  Transitioi 


Figure  6- 


PR0FIL07 
MSES  3.0 


(b)  Transitioi 


Figure  6.- 


PR0FIL07 


(b)  Transition  fixed. 


Figure  7.-  Concluded. 


0 


0  oc  =  -4.05 


(a)  a  =  -4.05°,  -3.04°,  -2.02°,  -1.00°,  and  0.02°. 

Figure  8.- Experimental  pressure  distributions  for  R=1.00xl06  and  M  =  0.11  with  transi¬ 
tion  free. 
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Figure  8.-  Continued. 
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(d)  a  =  11.20°,  12.20°,  13.20°,  14.19°,  and  15.18°. 


Figure  8.-  Continued. 
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(a)  R  =  0.50  x  106  and  M  =  0.05. 


Figure  9.-  Experimental  section  characteristics  with  transition  free,  transition  fixed,  and  rough. 
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(c)  R  =  1.00  X  106  and  M  =  0.11. 


Figure  9.-  Continued. 
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(d)  R  =  1.50  x  106  and  M  =  0.17. 


Figure  9.-  Concluded. 
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Figure  11.-  Variation  of  experimental  maximum  lift  coefficient  with  Reynolds  number. 
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Figure  12.-  Variation  of  experimental  minimum  profile-drag  coefficient  with  Reynolds  num¬ 
ber. 
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-  Theory  (PR0FIL07) 


a  =  3.17°  relative  to  the  x-axis 


(b)  C/  =  0.67. 
Figure  13.-  Continued. 


-  Theory  (PR0FIL07) 

-  Theory  (MSES  3.0) 

0 - 0  Experiment 


oc  =  9.85°  relative  to  the  x-axis 


0.5  x/c 

(c)  C/  =  C/ max. 


Figure  13.-  Concluded. 
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Figure  14.-  Continued. 


Theory (PR0FIL07) 
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Figure  14.- 


Theory (PR0FIL07) 


(a)  R  =  0.50  x  106. 


Figure  15.-  Comparison  of  theoretical  and  experimental  section  characteristics  with  transition  fixed. 
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(b)  R  =  0.70  x  106. 


Figure  15.-  Continued. 
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(d)  R=  1.50  x  106. 
Figure  15.-  Concluded. 
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R  =  0.50  x  106  and  M  =  0.05  with  Transition  Free 


a,  deg 

c  / 

-4.053 

-0.0354 

-3.035 

.0507 

-2.782 

.0625 

-2.528 

.0762 

-2.018 

.1279 

-.998 

.2206 

.021 

.3193 

1.038 

.4216 

2.055 

.5213 

3.073 

.6220 

4.090 

.7145 

5.107 

.8106 

6.123 

.8986 

7.139 

.9787 

8.154 

1.0565 

9.168 

1.1144 

10.179 

1.1610 

11.191 

1.2068 

12.198 

1.2225 

13.201 

1.2187 

14.197 

1.2014 

15.191 

1.1907 

16.178 

1.1197 

17.164 

1.0579 

cd 

cm 

0.011075 

-0.06872 

.009060 

-.06401 

.007873 

-.06255 

.007135 

-.05983 

.007448 

-.05801 

.008663 

-.05283 

.008899 

-.04943 

.009032 

-.04968 

.009232 

-.04951 

.009499 

-.04887 

.009974 

-.04718 

.010375 

-.04592 

.011043 

-.04440 

.011835 

-.04205 

.012971 

-.03864 

.014523 

-.03327 

.016671 

-.02843 

.019790 

-.02225 

.025501 

-.01644 

.020191 

-.01098 

.031392 

-.01255 

.051522 

-.01953 

.066080 

-.02171 

.081776 

-.02672 
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R  =  0.50x  106  and  M  =  0.05  with  Transition  Fixed 


a,  deg 

c  / 

-4.050 

-0.0232 

-3.033 

.0623 

-2.015 

.1631 

-.997 

.2640 

.021 

.3648 

1.038 

.4500 

2.056 

.5392 

3.071 

.6168 

4.088 

.7022 

5.104 

.7845 

6.119 

.8571 

7.135 

.9294 

8.149 

.9947 

9.162 

1.0393 

10.173 

1.0924 

11.183 

1.1351 

12.187 

1.1618 

13.189 

1.1795 

14.183 

1.1656 

15.178 

1.1568 

16.178 

1.1331 

17.167 

1.0888 

Cd 

cm 

0.010911 

-0.06797 

.009484 

-.06438 

.009492 

-.06287 

.009497 

-.06173 

.009923 

-.06093 

.010445 

-.05680 

.011342 

-.05353 

.012257 

-.04986 

.013251 

-.04644 

.014147 

-.04381 

.015149 

-.04019 

.016691 

-.03508 

.017447 

-.03106 

.020168 

-.02360 

.023030 

-.02016 

.027513 

-.01742 

.030099 

-.01719 

.030230 

-.01924 

.046598 

-.02469 

.065548 

-.03065 

.072738 

-.02562 

.089178 

-.03031 
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R  =  0.50  x  106  and  M  =  0.05  Rough 


a,  deg 

c  / 

-4.052 

-0.0280 

-3.033 

.0677 

-2.016 

.1533 

-.998 

.2325 

.022 

.3276 

1.039 

.4284 

2.057 

.5299 

3.074 

.6226 

4.091 

.7069 

5.107 

.7846 

6.121 

.8536 

7.136 

.9172 

8.149 

.9756 

9.161 

1.0265 

10.169 

1.0557 

11.174 

1.0687 

12.168 

1.0592 

13.153 

.9813 

14.141 

.9364 

15.130 

.9004 

16.120 

.8745 

cd 

cm 

0.011222 

-0.06842 

.010347 

-.06527 

.007926 

-.06189 

.007550 

-.05460 

.008682 

-.04976 

.009096 

-.04940 

.008861 

-.04924 

.009825 

-.04732 

.011361 

-.04389 

.012018 

-.03981 

.013582 

-.03591 

.015791 

-.03109 

.018609 

-.02658 

.023632 

-.02199 

.035178 

-.01738 

.050283 

-.01500 

.097468 

-.01949 

.130726 

-.02106 

.163616 

-.02766 

.196137 

-.03408 

.228290 

-.04114 
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R  =  0.70  x  106  and  M  =  0.07  with  Transition  Free 


a,  deg 

c  / 

-4.052 

-0.0313 

-3.035 

.0579 

-2.528 

.0914 

-2.274 

.1013 

-2.018 

.1298 

-.999 

.2308 

.019 

.3354 

1.037 

.4394 

2.055 

.5402 

3.073 

.6414 

4.091 

.7384 

5.108 

.8348 

6.124 

.9230 

7.141 

1.0098 

8.157 

1.0808 

9.170 

1.1404 

10.183 

1.1927 

11.194 

1.2218 

12.201 

1.2325 

13.201 

1.2209 

14.197 

1.2153 

15.186 

1.1512 

16.175 

1.1087 

17.159 

1.0284 

Cd 

cm 

0.010199 

-0.06762 

.009036 

-.06512 

.007172 

-.06308 

.006370 

-.06028 

.006800 

-.05926 

.007303 

-.05542 

.007602 

-.05586 

.007774 

-.05556 

.008035 

-.05449 

.008485 

-.05317 

.009025 

-.05178 

.009475 

-.05099 

.010352 

-.04943 

.011474 

-.04648 

.012620 

-.04107 

.014536 

-.03617 

.016573 

-.03046 

.020362 

-.02290 

.025432 

-.01572 

.029810 

-.01114 

.032923 

-.01719 

.041441 

-.01694 

.061829 

-.02329 

.076866 

-.02799 
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R  =  0.70  x  106  and  M  =  0.07  with  Transition  Fixed 


a,  deg 

c  / 

-4.050 

-0.0259 

-3.034 

.0618 

-2.015 

.1664 

-.997 

.2669 

.021 

.3633 

1.039 

.4586 

2.056 

.5382 

3.073 

.6278 

4.089 

.7077 

5.105 

.7917 

6.121 

.8715 

7.136 

.9436 

8.151 

1.0156 

9.165 

1.0771 

10.177 

1.1283 

11.187 

1.1659 

12.194 

1.1959 

13.195 

1.2107 

14.190 

1.2004 

15.183 

1.1564 

16.175 

1.1154 

17.161 

1.0535 

Cd 


0.009989 

-0.06635 

.009109 

-.06443 

.009002 

-.06364 

.009013 

-.06215 

.009027 

-.06057 

.009420 

-.05797 

.010008 

-.05323 

.011256 

-.05059 

.011873 

-.04731 

.012840 

-.04404 

.014378 

-.04077 

.015505 

-.03696 

.017166 

-.03344 

.019672 

-.02909 

.022347 

-.02374 

.028650 

-.01859 

.037090 

-.01710 

.028799 

-.01848 

.059526 

-.02327 

.054524 

-.02328 

.071827 

-.02607 

.086878 

-.03129 
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R  =  0.70  x  106  and  M  =  0.07  Rough 


a,  deg 

c  / 

-4.050 

-0.0256 

-3.033 

.0703 

-2.015 

.1626 

-.998 

.2455 

.019 

.3153 

1.036 

.3989 

2.053 

.4846 

3.070 

.5667 

4.087 

.6580 

5.103 

.7429 

6.119 

.8249 

7.134 

.8915 

8.147 

.9524 

9.158 

.9982 

10.164 

1.0168 

11.164 

1.0203 

12.155 

.9899 

13.143 

.9362 

14.132 

.8968 

15.117 

.8526 

16.107 

.8307 

17.096 

.8064 

Cd 

cm 

0.010351 

-0.06638 

.009815 

-.06548 

.009431 

-.06252 

.009210 

-.05832 

.009089 

-.05129 

.011085 

-.04659 

.012068 

-.04315 

.012488 

-.03908 

.013550 

-.03730 

.014803 

-.03560 

.016219 

-.03291 

.018491 

-.02764 

.022355 

-.02393 

.030426 

-.02012 

.048506 

-.01621 

.076856 

-.01795 

.126975 

-.02118 

.159901 

-.02439 

.192458 

-.03066 

.224645 

-.03994 

.256462 

-.04947 

.287912 

-.05855 
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R=1.00xl06  and  M  =  0.11  with  Transition  Free 


a,  deg 

c  / 

-4.054 

-0.0428 

-3.035 

.0541 

-2.526 

.1029 

-2.273 

.1161 

-2.019 

.1292 

-1.000 

.2351 

.019 

.3444 

1.038 

.4428 

2.056 

.5482 

3.074 

.6517 

4.092 

.7532 

5.109 

.8494 

6.126 

.9388 

7.143 

1.0208 

8.158 

1.0914 

9.172 

1.1548 

10.186 

1.2009 

11.196 

1.2339 

12.202 

1.2373 

13.201 

1.2532 

14.189 

1.1810 

15.179 

1.1267 

16.171 

1.0909 

17.163 

1.0630 

cd 

cm 

0.009723 

-0.06745 

.008415 

-.06475 

.007843 

-.06393 

.006670 

-.06224 

.005721 

-.05974 

.006201 

-.05838 

.006560 

-.05814 

.006794 

-.05530 

.007191 

-.05501 

.007569 

-.05503 

.008190 

-.05470 

.008829 

-.05355 

.009576 

-.05119 

.010720 

-.04678 

.012327 

-.04248 

.014316 

-.03766 

.016976 

-.02986 

.021070 

-.02396 

.026975 

-.01634 

.027641 

-.02104 

.034638 

-.02141 

.048480 

-.02294 

.063388 

-.02609 

.083186 

-.03184 
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R=1.00xl06  and  M  =  0.11  with  Transition  Fixed 


a,  deg 

c  / 

-4.053 

-0.0372 

-3.034 

.0588 

-2.015 

.1619 

-.996 

.2582 

.021 

.3433 

1.039 

.4268 

2.056 

.5187 

3.074 

.6120 

4.091 

.6970 

5.107 

.7826 

6.124 

.8669 

7.139 

.9459 

8.154 

1.0186 

9.168 

1.0826 

10.181 

1.1385 

11.191 

1.1798 

12.197 

1.1950 

13.189 

1.1640 

14.183 

1.1495 

15.178 

1.1318 

16.165 

1.0802 

Cd 

cm 

0.009609 

-0.06722 

.008415 

-.06452 

.008082 

-.06296 

.008760 

-.05906 

.009450 

-.05467 

.010473 

-.05047 

.011675 

-.04765 

.012520 

-.04519 

.013428 

-.04243 

.014088 

-.03963 

.015236 

-.03680 

.016778 

-.03355 

.018555 

-.02991 

.021246 

-.02599 

.025341 

-.02226 

.030733 

-.01815 

.035860 

-.01335 

.056752 

-.01554 

.034823 

-.02122 

.056535 

-.02615 

.072582 

-.03198 
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R=1.00xl06  and  M  =  0.11  Rough 


a,  deg 

c  / 

-4.054 

-0.0803 

-3.036 

.0090 

-2.018 

.1013 

-1.000 

.1975 

.018 

.2929 

1.036 

.3888 

2.054 

.4861 

3.071 

.5776 

4.088 

.6693 

5.104 

.7560 

6.120 

.8355 

7.135 

.9055 

8.148 

.9648 

9.159 

1.0064 

10.165 

1.0272 

11.158 

1.0225 

12.150 

.9789 

13.138 

.9244 

14.125 

.8718 

15.110 

.8266 

16.096 

.7968 

cd 

cm 

0.012093 

-0.05881 

.012189 

-.05437 

.011980 

-.05138 

.011892 

-.04874 

.012033 

-.04673 

.012225 

-.04488 

.012482 

-.04316 

.012490 

-.04143 

.013229 

-.03941 

.014070 

-.03755 

.015913 

-.03363 

.018444 

-.02941 

.023220 

-.02601 

.032627 

-.02089 

.051105 

-.01815 

.100854 

-.02530 

.134058 

-.02575 

.166891 

-.02892 

.199354 

-.03551 

.231445 

-.04507 

.263164 

-.05779 
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R=  1.50  x  106  and  M  =  0.17  with  Transition  Free 


a,  deg 

c  / 

-4.056 

-0.0549 

-3.037 

.0472 

-2.528 

.0952 

-2.019 

.1431 

-1.765 

.1539 

-1.000 

.2406 

.020 

.3440 

1.038 

.4533 

2.057 

.5613 

3.074 

.6615 

4.093 

.7646 

5.111 

.8600 

6.130 

.9549 

7.146 

1.0402 

8.162 

1.1161 

9.178 

1.1842 

10.191 

1.2380 

11.202 

1.2745 

12.209 

1.2921 

13.201 

1.2473 

14.189 

1.1781 

15.179 

1.1301 

16.170 

1.0938 

Cd 


0.008719 

-0.06723 

.008089 

-.06525 

.007919 

-.06459 

.006813 

-.06242 

.005245 

-.06059 

.005743 

-.05985 

.006021 

-.05725 

.006145 

-.05838 

.006496 

-.05844 

.007044 

-.05753 

.007477 

-.05641 

.008237 

-.05384 

.009397 

-.05144 

.010774 

-.04877 

.012585 

-.04489 

.015004 

-.03888 

.018020 

-.03340 

.022488 

-.02710 

.033325 

-.02302 

.025762 

-.02272 

.039779 

-.02372 

.053949 

-.02582 

.073737 

-.03108 
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R  =  1.50  x  106  and  M  =  0.17  with  Transition  Fixed 


a,  deg 

c  / 

-4.054 

-0.0472 

-3.036 

.0261 

-2.018 

.1137 

-.999 

.2087 

.019 

.3044 

1.038 

.4064 

2.056 

.5039 

3.074 

.6019 

4.091 

.6964 

5.109 

.7887 

6.126 

.8751 

7.142 

.9591 

8.158 

1.0358 

9.172 

1.1035 

10.185 

1.1612 

11.195 

1.2062 

12.196 

1.2093 

13.193 

1.2207 

14.175 

1.0919 

15.169 

1.0806 

16.156 

1.0269 

cd 

cm 

0.008935 

-0.06657 

.009869 

-.05906 

.010248 

-.05376 

.010621 

-.05031 

.010967 

-.04851 

.011310 

-.04688 

.011832 

-.04539 

.012369 

-.04395 

.013156 

-.04211 

.013990 

-.04008 

.015115 

-.03747 

.016472 

-.03476 

.018516 

-.03154 

.020699 

-.02808 

.024249 

-.02431 

.028151 

-.02062 

.044681 

-.02152 

.042494 

-.02840 

.042236 

-.02158 

.062191 

-.02844 

.081391 

-.03520 
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R=1.49xl06  and  M  =  0.17  Rough 


a,  deg 

c  / 

-4.057 

-0.1048 

-3.038 

-.0063 

-2.019 

.0988 

-1.000 

.1999 

.019 

.3007 

1.037 

.3990 

2.056 

.4973 

3.073 

.5907 

4.091 

.6837 

5.108 

.7709 

6.124 

.8528 

7.140 

.9268 

8.153 

.9912 

9.164 

1.0391 

10.170 

1.0636 

11.161 

1.0471 

12.150 

.9880 

13.138 

.9315 

14.123 

.8734 

15.107 

.8235 

16.092 

.7946 

Cd 

cm 

0.012571 

-0.05603 

.012085 

-.05395 

.011620 

-.05142 

.011537 

-.04942 

.011592 

-.04742 

.011992 

-.04593 

.012375 

-.04413 

.013060 

-.04209 

.014023 

-.03976 

.015361 

-.03692 

.017021 

-.03388 

.019846 

-.03044 

.024442 

-.02713 

.034944 

-.02397 

.056411 

-.02224 

.112803 

-.03053 

.145834 

-.03115 

.178490 

-.03326 

.210768 

-.03975 

.242671 

-.05106 

.274198 

-.06547 
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